A search for J PC = 1 + exotic mesons in the tt tt tt + and tt tt tt systems 
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A partial wave analysis (PWA) of the tt~7v~tt + and 7r _ 7r°7r° systems produced in the reaction 
n~p — > (3n)~p at 18 GeV/c was carried out using an isobar model assumption. This analysis 
is based on 3.0M tt~ty°-k events and 2.6M ir~ir~iv + events and shows production of the a2(1320), 
7T2(1670) and 7r(1800) mesons. An earlier analysis of 250K tv~tt~tt + events from the same experiment 
showed possible evidence for a J PC — 1~ + exotic meson with a mass of 1.6 GeV/c 2 decaying into 
pn. In this analysis of a higher statistics sample of the (37r)~ system in two charged modes we find 
no evidence of an exotic meson. 
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Quantum chrome-dynamics (QCD) predicts a spectrum 
of hybrid mesons beyond the qq bound states of the con- 
ventional quark model in which the gluons binding the 
quarks manifest their degrees of freedom. Lattice QCD 
and QCD-inspired models predict that the gluonic field 
within the meson forms a flux tube for quark separa- 
tion greater than 1 fermi yij. When the flux tube is in 
its ground state conventional mesons emerge. When the 
flux-tube is in its excited states hybrid mesons emerge. 
For conventional mesons the spin (J), parity (P) and 
charge conjugation (C) quantum numbers of the qq sys- 
tem are those of a fermion-antifermion system: J = L+S 
where L is the angular momentum between the quarks 
and S is the total quark spin; P = (— 1) +1 ; and 
C = (-l) i+s . Thus J PC = Q+-,l-+,2+-,... are 
not allowed. These exotic quantum numbers are allowed 
when the quantum numbers of the excited flux-tube are 
included and are thus a signature for hybrid mesons. The 
fundamental interest in the spectrum and other proper- 
ties (decay modes and widths) of hybrid mesons is that 
flux tubes are thought to be responsible for the confine- 
ment of quarks and gluons in QCD. 

In 1998, the E852 collaboration reported evidence for 
the 7Ti(1600), a J PC — 1 h exotic hybrid meson with a 
mass of 1.6 GeV/c 2 decaying into pir. The original publi- 
cation 2] of this result was followed by another |3| giving 
more details of the analysis based on 250,000 events of 
the reaction n~p — > ir~ir~ir + p at 18 GeV/c collected 
in the E852 experiment at Brookhaven Lab in 1994. A 
search for this state by the VES collaboration in the reac- 
tion tt~A — > 7r~7r~7r + A interactions at 37 GeV/c jj] did 
not confirm this state. The partial wave analysis (PWA) 
procedure used to extract these results was based on the 
isobar model assumption whereby the final state 3tt sys- 
tem is reached via an intermediate state of a di-pion res- 
onance (e.g. a, / (980), p(770), / 2 (1275), p 3 (1690)) and 
a bachelor n. 

This paper reports on the analysis of additional data 
collected in 1995 including 3.0M events of the reac- 
tion ir~p — > 7r~7r 7r°p and 2.6M events of the reaction 
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FIG. 1: Acceptance-uncorrected (a) 7r _ 7r°7r° and (b) 
■k~tt~tt + effective mass distributions along with acceptance 
functions and the (c) 7r~7r°7r° and (d) tt~it~tt + Dalitz plots 
for events in the 7r2(1670) mass region. 



n~p — > ir~ir~ir + p. We compare the isobar model PWA 
of these two large data sets with each other and with 
the analysis presented in 0, • The comparison of the 
two (37r)~ modes provides powerful cross checks. Any 
resonance decaying to (p7r)~ should decay equally to 
p~7r° — ► (7r~7r°)7T and p°ir~ — ► (7r + 7r - )7r~ and thus ap- 
pear with equal probabilities in the two modes. Similarly, 
any resonance decay to /27r~ ought to appear twice as of- 
ten in (7r + 7T~)7r~ as (7r°7r°)7r~. And these same isospin 
rules dictate relationships between the two modes for all 
other decay sequences. Since the ir~ir~ir + and 7r~7r°7r° 
modes rely on different elements of the detector, any mis- 
understandings in the acceptance would affect the two 
modes differently and inconsistencies would result. As 
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FIG. 2: Comparison of the (a) H(201) and (b) H(420) mo- 
ments as computed directly from data and from PWA fits for 
the low and high wave sets for the ty~-k°tt channel. In (c) 
the A 2 (differences squared divided by errors squared summed 
over all mass bins and divided by the number of mass bins) 
for various moments. 



FIG. 3: Comparison of the (a) H(201) and (b) 11(420) mo- 
ments as computed directly from data and from PWA fits for 
the low and high wave sets for the ir~ir~Tv + channel. In (c) 
the A 2 (differences squared divided by errors squared summed 
over all mass bins and divided by the number of mass bins) 
for various moments. 



the results will show, this is not the case. 

The data reported on in this paper were collected in 
the E852 experiment at the Alternating Gradient Syn- 
chrotron (AGS) at Brookhaven Laboratory (BNL) using 
the multi-particle spectrometer (MPS) augmented with 
a forward electromagnetic calorimeter (LGD) consisting 
of 3000 lead glass blocks. An 18.3 GeV/c ir~ beam was 
incident on a liquid hydrogen target. Details about the 
E852 experiment are given elsewhere 

The 1995 run of E852 collected approximately 146M 
(265M) triggers requiring one (three) forward going 
charged particles. For those triggers requiring only one 
forward track, it was also required that the trigger pro- 
cessor detected an energy deposition pattern in the 
LGD consistent with an effective mass greater than a 
single 7T°. The kinematic fitting program SQUAW 7] 
was used in identifying approximately 73M (79M) recon- 
structed events that were consistent with the topology 
7r~p — > (3ir)~p. This constrained fit used the measured 
track and beam momenta, the measurements from the 
LGD and constraints on the recoil particle mass (= m p ) 
and di-photon effective masses (two pairs = m^) to re- 
fine the measured final state particle 4-momenta. In the 
it~it it topology, 13. 8M events survived this selection 
while 16. 8M events of the ir~Tr~ir + topology survived. 
Further cuts on fit confidence level, vertex position and 



correlation between the 37r system and a recoil track 
produced the final data set analyzed here: 3,025,980 
7r~7r°7r°events and 2,585,852 7r~7r~7r + events. 

Figures [I] (a) and (b) show the acceptance uncorrected 
7r _ 7r°7r° and 7r - 7r~7r + effective mass distributions respec- 
tively. Enhancements in the <zi(1260), <22(1320) and 
7T2(1670) mass regions are observed. The acceptance in 
37T effective mass is also shown in Figures ^ (a) and (b). 
Dalitz plots for the 7r _ 7r°7r° and ir~ir~n + systems in the 
7r 2 (1670) mass region are shown in Figures [I] (c) and (d) 
respectively. These Dalitz plots clearly show rich struc- 
ture and bands corresponding to decays into 7T/o(770) and 
7r/2(1275), motivating application of the isobar model. 

The PWA software used in this analysis consisted of 
programs developed at BNL and Indiana University (IU) 
Q . The formalism is similar to that used in 0, Q but 
the details of handling the data and doing the the PWA 
fits differ. The IU software was optimized for running 
on a 200-processor computer cluster (AVIDD) 9] allow- 
ing systematic studies involving many PWA fits, varying 
isobar parameters and using different wave sets. Results 
from the two programs were compared with each other 
and found to be consistent. The IU programs were also 
used to analyze the data presented in |2],|3j] and produced 
consistent results. 

The analysis is performed in the Gottfried- Jackson 
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FIG. 4: PWA results for the ty it°w° channel as a function 
of 37r effective mass: (a) sum of all waves and the background 
wave; (b) 2++ wave; (c) 2~+ wave; (d) A$(2 ++ - 2 _+ ); (e) 
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For (b) through (f ) the results 



for the low wave and high wave sets are shown. 



FIG. 5: PWA results for the n tt n channel as a function 
of 37r effective mass: (a) sum of all waves and the background 
wave; (b) 2 ++ wave; (c) 2"+ wave; (d) A3>(2 ++ - 2" + ); (e) 
1~+ wave; (f) A$(2 ++ -l"+). For (b) through (f) the results 
for the low wave and high wave sets are shown. 



frame. Within the isobar formalism each partial wave 
is characterized by J PC M e , (SL), where J, P and C are 
the spin, parity and charge conjugation of the 3-7T sys- 
tem, respectively, M is the spin projection along the z 
axis and e represents symmetry of the 3tt system under 
reflection in the production plane. Finally S is the spin 
of the isobar and L is the relative orbital angular momen- 
tum between the isobar and the bachelor 7r. Except for 
the S = wave, all 7T7t resonances are parametrized using 
relativistic Breit-Wigner functions with Blatt-Weisskopf 
factors. 

Following 01 a sufficient set of partial waves was de- 
termined by sequentially removing waves from a parent 
set containing waves with J < 4, M < 1 and S < 3 and 
examining the resulting change in likelihood. If no sig- 
nificant change in likelihood was obtained, the wave was 
removed. The exotic 1 ^ partial waves were kept even 



though they would have been chosen for removal by the 
above criteria - the existence of signals in these waves 
is discussed below. Based on these criteria the compiled 
set includes 35 waves and a background wave- we refer 
to this as the high wave set. In comparison, the wave set 
used in the analysis described in |3| included 20 waves 
and a background wave - we refer to this as the low wave 
set. The 35 wave set included three additional 1 ++ , seven 
additional 2 ^ and seven additional J > 3 waves This 
high wave set did not include two negative reflectivity 
partial waves (1 ++ 1~,(10) and 2 ^1^,(20)) that were 
included in the low wave set PWA. 

Another arbiter of wave set sufficiency is the compar- 
ison of moments, H(LMN), of the D^ IN (fl) functions 
as calculated directly from the data and as calculated 
using the results of PWA fits. We define H(LMN) = 
J I(H,)D^ IN (il)dn where represents the Euler angles 
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of the 37r system [ll| and the intensity /(O) is deter- 
mined directly from experiment or computed using the 
results of the PWA fits. In FigureHwe show the H(201) 
and H(420) moment comparisons as a function of 3ir mass 
for the 7r - 7r 7r° channel for low wave and high wave set 
PWA fits. We also show the difference between data and 
PWA calculations of the moments as A 2 (summing dif- 
ferences squared divided by errors squared summed over 
all mass bins and divided by the number of mass bins) 
for various moments. Similar plots are shown in Figure|3| 
for the 7r - 7r~7r + channel. The moments calculated us- 
ing the PWA results from the high wave set have better 
agreement with experimental moments. But there are 
moments (such as H(201) - shown in part (a) of Fig- 
ure |3J) for which agreement is not achieved and this is 
most likely due to the inherent inadequacy of the isobar 
model in describing the underlying production mecha- 
nism. 

Figures 0] and [5] show the PWA fit results for the 
ir~ir°ir° and 7r~7r~7r + channels respectively as a function 
of the 3ir effective mass. In part(a) of each figure the 
sum of the waves is shown along with the background 
wave. The background wave is added incoherently with 
the other waves and is isotropic in all decay angles. The 
fitting procedure constrains the sum of all waves to de- 
scribe the observed 37r effective mass. Parts (b), (c) 
and (d) of each figure show the mass dependence of the 
2++1+, (12) and 2^+0+, (20) waves and their phase dif- 
ference. Results are shown for both the low wave and 
high wave set. To within 10% the acceptance corrected 
yields of in the two 37r modes are equal, consistent 
with expected decays into p~ir° and p°ir~ . And the yield 
for is higher by a factor of two for the /2(1275)7r~ 
modes when f% — > ir + ir~ compared to fa — ► 7r°7r°, con- 
sistent with expectations from isospin. The 2 ++ and 2 v 
waves and their phase difference are consistent with two 
interfering Breit-Wigner line shapes. For the 2 h wave 
in the 7r - 7r~7r + channel (Figure EJc)) we also show the 
published 2 h wave from the PWA of diffractive 3-7T pro- 
duction in tt~p interactions at 63 and 94 GeV/c from the 
ACCMOR collaboration [l^. 

The 1 h l + , (11) wave and its interference with the 
2 ++ wave are shown in parts (e) and (f) of each figure. 
When the low wave set is used in the fit, an enhance- 
ment is observed for the exotic 1 ^ wave in the 3it mass 
region around 1.6 GeV/c 2 , consistent with the observa- 
tion of 0, 0| _ *■ e. we obtain the same results when we 
use the same wave set. But the enhancement disappears 
when the high wave set is used. This effect is observed 
for both 37T channels. And the phase (see part (d)) of the 
exotic wave relative to the dominant 2 ++ wave is similar 
for both wave sets and dominated by the phase of the 
a2(1320). The amplitudes of the negative reflectivity ex- 
otic waves show the same behavior, disappearing as the 
fit improves. 



The masses and widths of the p and the /2(1275) were 
varied and the maximum likelihood was achieved using 
current best values for these quantities 0|. Different 
parameterizations of the irir S-wave were also examined 

El Q • While small changes in the likelihood were 
observed, no systematic changes in the PWA results oc- 
curred. 

The PWA results presented here were for data with 
0.18 < \t\ < 0.23 (GeV/c) 2 where t is the square of the 
momentum transferred from the incoming 7r~ to the out- 
going (37r)~ system. PWA fits were also carried out for 
12 other regions in t and the results are similar to those 
presented here. 

In summary, there is no evidence for an exotic J PC — 

1 h meson in the 37r system and in particular for the 

7Tl(1600). 



ACKNOWLEDGMENTS 

This work was supported by grants from the 
United States Department of Energy (DE-FG- 
91ER40661/Task D and DE-FG0287ER40365) and 
the National Science Foundation (EIA-01 16050) for 
AVIDD. 



* Also with the Nuclear Theory Center, Indiana University 
[1] K. J. Juge, J. Kuti and C. Morningstar, hep-lat/0401032 
(2004). 

[2] G. S. Adams et. al. Phys. Rev. Lett. 81, 5760 (1998). 

[3] S. U. Chung et. a l Phys. R ev. D65, 072001 (2002). 

[4] I. Kachaev et. al, |hep-ex/0111067| (2001). 

[5] T. Adams et. al, Nucl. Inst. & Meth. A368, 617 (1996); 
Z. Bar-yam et. al, Nucl. Inst. & Meth. A386 235 (1997); 
and S. Teige et. al, Phys. Rev. D59, 012001 (1999) 

[6] B. B. Brabson et. al, Nucl. Inst. & Meth. A332 419 
(1993); R. R. Crittenden et. al, Nucl. Inst. & Meth. 
A387, 377 (1997); J. Gunter, Ph. D. Thesis, Indiana 
University (1997); and R. Lindenbusch, Ph. D. Thesis, 
Indiana University (1997). 

[7] O. Dahl et. al, SQUAW kinematic fitting program, 
Group A programming note P-126, U. California Berke- 
ley (1968). 

[8] Details about the IU software and PWA can be found at 

the website: dustbunny.physics.indiana.edu/3p_paper. 
[9] AVIDD (Analysis and Visualization of Instrument Driven 

Data) Website: kb.indiana.edu/data/almb.ose.help. 
[10] C. Daum et. al. (ACCMOR Collaboration), Nucl. Phys, 

B182,269 (1981). 
[11] G. Ascoli and H. W. Wyld, Phys. Rev. D12, 43 (1975). 
[12] S. Eidelman et. al, Phys. Lett. B592, 1 (2004). 
[13] K. L. Au, D. Morgan and M. R. Pennington, Phys. Rev. 

D35, 1633 (1987). 
[14] J. A. Oiler and E. Oset, Phys, Rev. D60, 074023 (1999). 



